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1. Introduction 

Integrable hierarchies of nonlinear partial differential equations (PDE) have been 
vigorously studied from the perspective of physics as well as mathematics. Among 
them the Kadomtsev-Petviashvili (KP) hierarchy and its variants appear in many areas 
of theoretical physics. In particular, the supersymmetric extensions of the KP hierarchy 
play important roles in non-perturbative superstring theories and the connections 
are suggested between the dispersionless limit of the KP hierarchy and topological field 
theory g i. 

Concerning with the construction of the KP hierarchy in the Lax formalism, the 
non-commutative algebra of pseudo- or micro-differential operators play fundamental 
role IQ. For the standard KP hierarchy, the associated pseudo-differential operator can 
be regarded as an ordinary integral operator, which enjoys the generalized Leibniz rule. 
The aim of this paper is to inquire into the practicability of the extensions of the KP 
hierarchy by introducing fractional order pseudo-differential operators. In this respect, 
we recall that the survey of the fractional order integration and differentiation is known 
as fractional calculus. 

The fractional calculus, which usually stands for the differentiation and integration 
of arbitrary order so the terminology is somewhat misleading, has a long and rich 
history [^, The standard definition of the arbitrary order integration/differentiation 
is mostly given by the so called Riemann-Liouville integral these days. Although 
the fractional calculus has been studied well in mathematics, it is not an ordinary 
mathematical tool in the theory of integrable systems at present. Apart from integrable 
systems, there are many applications of fractional calculus in physics; for example, one 
of the present authors analyzed the supersymmetric field theories through half-order 
differential operators [^], other important works in the subject are performed on non- 
differential evolution equations, chaotic dynamical systems, material physics, and so on 

i- 

In the present paper we consider extensions of the KP hierarchy by introducing the 
fractional order integral/differential operators as pseudo-differential operators, which 
should be interpreted as the power roots of ordinary integration/differentiation; the 
situation is similar to the supersymmetric extensions of the KP hierarchy |10|, [l^. 



12|, where the square-root of integral/differential operators are brought in through 
superspace formulation. In contrast, we extend the KP hierarchy by making use 
of fractional order integral/differential operators with respect to purely "bosonic" 
variables, for which the relevant non-commutative algebra is the generalized Leibniz 
rule of fractional order. We will see in the following that the extension of the KP 
hierarchy by half-order integral/differential operators leads to a hierarchy being similar 
to that of supersymmetric extension, as expected. 

This paper is organized as follows. In the next section we give a very brief review 
of the Lax formulation of the KP hierarchy and its supersymmetric extension for the 
purpose of determining notation. In section 3 we make an attempt to generalize the KP 
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hierarchy by fractional order integral/differential operators, and find the formulation 
works consistently. In section 4 we introduce the Riemann-Liouville integrals as a 
candidate for the pseudo-differential operators of fractional order, which supply the 
generalized Leibniz rule being used in section 3. The final section is devoted to 
concluding remarks. 



2. The LcLx formulation of the KP hierarchy 

In this section we give a sketch of the Lax formulation of the standard KP hierarchy, its 
fc-reduction and supersymmetric extensions, to fix the notation throughout the present 
paper. 



2.1. The standard KP hierarchy 

The KP hierarchy within the framework of Lax formulation is generated by the 
non-commutative algebra of the pseudo-differential operator with respect to an 
independent variable x, which acts on a function through the generalized Leibniz rule. 



d~'-f = Yl[^)f^'^d~^-'- (1) 

Here we consider the case of integer j, the order of pseudo-derivative or integral, although 
the formula (P is valid for non- integer j. We define the Lax operator of the (one- 
component) KP hierarchy by, 

oo 

Ckp = d +Y,^3+id-\ (2) 

i=i 

where m/s are dependent variables of space x and time variables being introduced below. 
The coefficient of d*^ can be set zero without loss of generality. We assign the degree of 
the differential operator d one, standing for deg[9] = 1, and assume that all the terms in 
the Lax operator (||) have equal degree, i.e., deg[uj\ = j. This assignment of the degree 
is naturally justified by the tau-function formalism. Introducing infinite directions of 
"time" t = (^1,^2, ^3, • • •) with deg[t„] = —n, we may consider the Lax equations, 

fir 

^^ = [B^Xkp] (n = 1,2,3,...). (3) 

If we define the n-th "Hamiltonian" Bn by the non-negative power part of d in the n-th 
product of the Lax operator (g), denoting Bn '■= {C^p)+ = (£^p)>o, we will obtain an 
infinite tower of nonlinear PDE's, the standard KP hierarchy. Note that the lowest time 
variable ti should be identified with the space variable x due to the first Lax equation. 
The lowest PDE, the KP equation, is obtained by comparing the coefficients of ^"■''s in 
each side of (§) for ^2 and t^ developments of U2, and and eliminating the ^3 and 
M4, 



Ady^ 



du 1 ; 

-u — 6uu 

dt 4 



(4) 
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where u := ^2, U '■= t2 and t := t^, and the prime is the derivative with respect to 
x{= ti). 

2.2. The k -reduction 

The KP hierarchy is an unconstrained system in the sense that the dependent variables 
Uj are independent of each other in the Lax operator @. This independence is not 
necessary: we can impose constraints between dependent variables without loss of 
consistency. The most familiar is the /c-reductions of the KP hierarchy for an integer 
k > 2, for which the constraints are -C^p — ^k, ^-c, all the coefficients of negative 
powers in d of C^p are zero, 

{^KP).^ = 0, (5) 

where m = 1, 2, . . .. This is equivalent to the tik independence of the system, for a natural 
number /. For example, the 2- and 3- reduction lead to the KdV and the Boussinesq 
hierarchy, respectively. For later purpose, we make a trivial remark that the reduction 
conditions (^) are compatible with the Lax equation: the conditions are invariant with 
respect to the time evolutions, 

ark 



dtn 



= 0, (6) 
because the Hamiltonians have only derivatives, i.e., non- negative power terms in d. 

2.3. Supersymmetric extensions 

Supersymmetric extensions of the KP hierarchy (SKP) are vigorously studied by both 
mathematicians and physicists. In particular, they appear in the context of superstring 
and/or quantum gravity theories The first supersymmetric extension was done 
by Manin and Radul [§, referred to MRKP^, in which the differential operator in 
superspace, i.e., the superderivative, and its inverse, 

play the parallel role of d and in the standard "bosonic" KP. Here, ^ is a Grassmann 
odd variable, accordingly the square of D turns out to be ordinary derivative, 

in other words D can be regarded as a square root of d. According to the superspace 
formalism, superfields $j play the role of dependent variables in the MRKP, whose Lax 

% Another formulation of SKP is given in |ll|. 
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operator of the MRKP is defined as, 



(9) 



Besides the bosonic time variables t, infinite fermionic time variables (ri,r2, . . .) must 
be introduced. Consequently we observe that both the bosonic and fermionic time flows 
of the superfields make up a system of super-differential equations. 

To make a comparison with another extension of the KP hierarchy considered in 
the following section, we exhibit the lowest degree bosonic time flows of the MRKP, 
which is given by the Lax equation of even order, 

dCuR 



dtr 



MR 



where the Hamiltonian is the standard one: 



2n 



2n \ 
MR) 



(10) 



In addition, there exist 



fermionic time flows given by odd order Lax equation certainly, we do not need them. 

One can show the lowest 
[1), which can be given in 



however, in the present consideration, for the detail see ^ |T 
degree equation of (jT^) is an extension of the KP equation 
the component form 

Sd'^u \du 1 ,„ „ , 3 
'di 



4dy2 



dv 
di 



4 



r 



— 3uu' 



-v"v 



(llfl) 

(116) 
V + 6u, and 



where the bosonic variable u and the fermionic one v are defined by -D$2 
t := ts and y := ^2- 

Besides the MRKP, various types of supersymmetric extension of the KP hierarchy 



are considered [0, |T5[ . For example, a non-standard Lax equation by Brunelli and Das 
leads to an extension of the KP equation of the following form. 



3d^u 



Ady^ 



du 
'di 

dt 4 



3uu' 



-V V 



-V 



^ dv 3 dv 
dx — —!!■ 



dy 



2 dy^ 



^[nv) - -u 



^ dv ^ 3 , rdu , 
—dx + -V / —dx 
dy 2 J dy ^ 



'12a) 



'12b) 



where, similarly to the MRKP, u and v are bosonic and fermionic variables, respectively. 
In contrast to (|lla|) and ( |11&|) , there appear non-local terms in these coupled equations. 



3. Extensions of the KP hierarchy by fractional order integral operators 

This section provides the extensions of the standard KP hierarchy by fractional order 
integral operators, which is the main topic of the present paper. 

3.1. Extended Lax operator 

Recall that the Leibniz rule (|lD is applicable when the order j of "integral" is an arbitrary 
real (or complex) number. It will be interesting to consider the case when the Lax 
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operator includes fractional order integrals, and then, to inquire whether the system 
gives a consistent hierarchy or not In the following consideration, we accept the 
axiom that the fractional order integral operators exist and also its exponential law 
Q-'^Q-o = Q-{i-+j) holds for fractional i and j, for a while. We will see the Riemann- 
Liouville integral of fractional order enjoys these requirements in the next section. 

3.1.1. Extension by the half-order integrals For the simplest case of an extension of the 
KP hierarchy, we consider the Lax operator including the half-order integrals in addition 
to the Lax operator (H). We restrict ourselves to the case that the highest order term 
is d as in the KP. Accordingly, we define the most general half-order integral operator, 

where fm's are the dependent variables of degree m/2. We have set the "differentiation" 
term d^^"^ to be absent: this resulted from the Lax equation defined below. We remark 
that the Lax operator composed only of the half-order integrals (0) itself does not 
produce any consistent hierarchy, because its products does not close in the half-order 
integral operators: we need integer order integral/differential operators to close the 
algebra. With this definition, we consider the following Lax operator, 

Cl/2 = CkP + Mi,2, (14) 

and the standard Lax equation for the flows with respect to the time t = (ti, ^2, • • 

^ = [i3„,A/2]. (15) 

If we take the standard definition of the Hamiltonian, Bn '■= (>C"y2)+5 whose lower degree 
sequence is, 

Bi = d (16a) 
B2 = d^ + 2v3 d^l'^ + 2u2 (166) 
^3 = 5^ + 3t;3 (9^/2 + 3^29 + 3(t;5 + v^)d'l'^ + 3^3 + 3^2 + St^g, (16c) 

then we find closed coupled nonlinear PDE's, an extended KP hierarchy by the half-order 
integrals, hereafter eKPi/2- Other "non-standard" definitions of Bn such as ('C"y2)>i/2 
cause inconsistency. To show the consistency of the system, we derive the lowest degree 
coupled PDE from (|1^), the extended KP equation by the half-order integral, i.e., the 
eKPi/2 equation. Just like the original KP equation (^, we need the first two non-trivial 
equations of (P!3|). Each coefficient of the negative powers in d of 



dH\i2 
dt2 

is, 



[B2,Cy2], (17) 



d-'/' : ^ = 2v', + v'^ (18a) 

Here we restrict ourselves to rational j; if j is irrational, the Lax equation could not give a closed 
system, see the following argument. 
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du2 



d-' : 



dy 



+ U2 + 2t>3t>3 



whereas of 



IS, 



d- 



dCi/2 
dt. 



du3 
dy 



9-1/2 : ^ = 3t;; + 3< + < + 6(W 



d 



-1 



du2 
'dt 



3^4 + 3^3 + u'2 + 6U2U2 + 6(f3f5)' 



+ 2(^3 +^3^3)- 
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(186) 
(18c) 
(18d) 

(19) 

(20 a) 

(206) 



Eliminating the dependent variables u^, u^^ and W7, we find the coupled PDE with 
non-local term, 



3d^v 



dt A 8^ ^ 4 y 



dy 4 dy 



where u := U2 and v 



-3(H 



(21a) 



(216) 



^3. As expected, (|21q|) reduces to the KP equation (0) 
when V is absent. We observe the resemblance between ( pi 4 ), (pl6|) and the MRKP 
equations ( p.lc|) , (|116|) or the non-standard SKP equations (|12q|) , (|12£|) , however they 
are not exactly identical. This resemblance obviously comes from the fact that the 
derivative in superspace can be read as a square root of the derivative, which fact is 
formally equivalent to the feature of the half order derivative d^^"^. In contrast to the 
supersymmetric extensions, the extension considered in this section works without using 
Grassmann numbers. 



3.1.2. Extension by the 1/N-th order integrals Having observed the extension by the 
half-order integrals is successful, we now consider more generic extensions by the 1/A^-th 
order integrals (A^ = 3, 4, . . .), eKPi/7v hierarchies. In these cases, we need to introduce 
integral operators d~^^^ , d'"^^^ , . . . , d^^^^^^^^ simultaneously to give a consistent Lax 
equation, since we have to close the commutator algebra in the Lax equations under 
the axiom d~^d~^ = d~^~K For N = p, a prime number, there appears a new system 
coupled to the KP hierarchy. For example, we give an outline of the = 3 case, in 
which the Lax operator should be made up of, 

Ci/3 = CKP + M1/3 + M2/3, (22) 
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M,/3 = + d-^l^ + w^o d-"^ + • • • , (23) 

M2/Z = d-''/^ + ws d-''/^ + wri d-^/^ + • • • , (24) 

and deg[-u;TO] = m/3. We observe that the standard Lax equation and the definition of 
the Hamiltonian similar to the former case give a consistent hierarchy of coupled PDE's. 
One can see the lowest coupled PDE arises from the first two non-trivial Lax equations. 
Each coefficient of d in 

^^^/^^[i3„A/3], (25) 



IS, 



whereas of 



IS, 



dy 

d-'/':^ = 2w', + w'i + {wiy (266) 

d-^ : -p-^ 2u'^ + 4' + 2(^4^5)' (26c) 

a-V3 . ^ ^ 2w[^ + < + {wiy + 2{w^U2)' (26d) 

9-^/^ : ^ = 2wn + w'l + ^^^7^ + ^Wiw'^ + 2^5-^2 - '^Wiw'i (26 e) 
"^Ms ^ , // r. / 10 / 4 ,8 , 

d : — — = 2^4 + % + 2U2U2 + — ty8«^4 + T7^4«'8 + 77^t^7«^5 

a^/ 3 3 3 

2 , 4 ,, 2 ,, , ^ 

+ gWgWy - -W5W4 - -W4W5, (26/) 



<9£i 



/3 



dt. 



= [B3,A/3], (27) 



. ^ ^ + 3^^' + + 6(«;4ii2)' + ^{wiy (28a) 

. ^ _ ^ ^ ^/// ^ 6(W5«2)' + 6(W7W4)' 

at 

+ 2w;^ + 2W4W4 (286) 

d'^ ■ = 3^4 + 3^3 + U2 + 6U2U2 + 6(wsW4y + diw^w^y 

at 

+ w^w" + + 3wlw[. (28 c) 

These are nine equations for the nine dependent variables so that we can combine them 
into the coupled PDE of U2, W4 and w^. 

For N being a composite number, we observe that the new system is coupled to 
the system coming from the prime factors of N. For example, the eKPi/4 system is a 
new system coupling to the eKPi/2 system given above. 
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Finally we should remark that the introduction of pseudo-derivative of irrational 
order does not make a finite closed system: we need uncountable number of additional 
M's like and (0). 



3.2. The conservation laws 

Since the eKPi/Ar hierarchy is constructed within the framework of Lax formalism, we 
expect that the system is integrable a priori. In fact, we observe that there are infinite 



conservation laws, which can be derived by standard procedure ||2^ for the Lax operator 
under consideration, say then we find, 

ARes(£;) = p;_^, (29) 

where the residue is defined as (£ and Pm,n is a differential polynomial of Uj and Vj. 
One can see that the presence of the fractional order integral operators do not modify 
the formula (^91). Hence we expect the existence of many special solutions to the eKPi/7v 



hierarchies, just like the solitons in the original KP. 

In addition to these conserved charges with integer degree, we have another set 
which has non-integer degree. For concreteness, we consider the eKPi/2, in which there 
exist conserved charges with degree k + 1/2 {k = 0, 1, 2, . . .): we can find the charges 
come from, 

Res(/:^;/), (30) 
where the square root of the Lax operator is constructed by usual procedure, 

^i/l = + l^sd-' + l^2d~l + ^iv,- \v',)d-' + ■■■. (31) 
Although we can construct the charges with half-integer degree, there does not exist a 

k+- 

consistent time flow generated by the Hamiltonian B|^_^_l := jC^^^ . 

In general, for the eKPi/7v (A^ > 3) we will find the existence of additional sequences 
of conserved charges. 



3.3. The k-reduction of the eKPi/^ 

In this subsection we make a comment on the fc-reduction of the cKPi/at, the truncation 
of the tik flow. Unlike the standard KP hierarchy, the reduction condition Cyj^ = Bk 
does not work in the eKPi/jv hierarchies due to the property of fractional integrals. For, 
the compatibility (|) between the reduction condition and the Lax equation does not 
hold when the fractional order "derivative" operators are present in the Hamiltonian Bk- 
This comes from the fact that the Leibniz rule for non-integer order is not a finite sum 
even if the order j is positive, hence the right hand side of the corresponding equation 
to (^ induces negative terms in d, i.e., 

i[Bk,Bk'])_^^0, (m = 1,2,3,...). (32) 
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For example, the commutator of the Hamiltonians B2 and for the eKPi/2 gives, 

{[B,M)-^^ = \{v,u',y -{vly (33) 

as well as the coefficient of 9"™/^ (m = 2,3, . . .). Although it is not clear at present 
whether one can impose consistent reduction conditions on the Lax operator of the 
eKPi/^, we may consider the truncated system by hand. 

Here we present a remarkable fact that there exists an algebraic solution to the 



2-reduction of the eKPi/2 equation ( |21o|) and (|21£| ), which should be referred eKdVi/2 
equation, 

^-^^=\u"' + ?,uu'-\{vy, (34a) 

1^ = \v"' + Kuvy. {Uh) 
Notice the strong resemblance between the supersymmetric KdV equation pO[ and 



( |34?j ) . The solution is an extension of the rational solution to the KdV equation: 

uixA) = — — - 777 + -, (35a) 

^ ' ' 16(x + ct)2 3' ^ ' 

, , , VT65 1 /o^;\ 

24 [x + ctyi^ 

where c is a constant with deg[c] = 2, which is required by keeping the correct degree 
of the dependent variables. 

4. The Riemann-Liouville integrals of fractional order 

So far we constructed the extensions of the KP hierarchy by introducing the pseudo- 
differential operator of fractional order, the fractional integral. However, we have treated 
such operators as only generators of a non-commutative algebra. In the construction 
of the standard KP hierarchy, we use the Leibniz rule of negative order derivatives (|I|), 
which can be realized by the iterative use of integration by parts, e.(?., for j = 1, 

d-\fg) = fg dx = fG^'^- f'G^'Ux 

= fG(^) - f'G(^) + f'G^'^dx 

= fQi^) _ f'G^^) + - ■ • • 

00 

= ^{-lff^^^G^''^^\ (36) 

fc=0 

where G^^^ is fc-times indefinite integral of g. The cases of higher order j can be shown 
by multiplicative operation of (pG]). Hence we accept the statement that the negative 
order derivative operator is equivalent to indefinite integral operator. 
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For the case of fractional j, how can we reahze the formula (|1|)? To find out the 
appropriate fractional integration/derivation on a function assumed to be existing in 
the last section, we recall the Riemann-Liouville integral of order a G C for 3?e a > 0, 

r/(a;) := -i- f^x - zr~'f{z)dz, (37) 
r(«) Ja 

where f{x) is assumed to be locally integrable and rapidly decreasing on the lower 
boundary a. We realize that when 3?e a < 0, a ^ —N, (^) should be read as, 

rf{x) := = — (38) 

where n is the first integer of n + 3?ea > 0. The assumption for f{x) guarantees the last 
equality in (^), which means that the Riemann-Liouville integral commutes with the 
ordinary derivative. Hence we understand that the definition ( p7D can be extended by 
analytic continuation to the whole a, for the detail see Appendix. One can also show, 
when a tends to an integer, the Riemann-Liouville integral turns out to be the ordinary 
integration/differentiation, as it should be. In addition, we will observe in Appendix 
that the Riemann-Liouville integral complies with the exponential or additive index law, 
_ /Q!+/3^ which is assumed in the derivation of the eKPi/^v hierarchy. 
We exhibit the action of the Riemann-Liouville integral on some functions, 

where A > and, 

ja M = ^il^+'^) M+« (40) 
r(a + /i + l) ' ^ ^ 

where the operand, x^, is defined zero if x < 0. 

In Appendix we observe that the definitions ( P7| ) and ( pS] ) lead to the Leibniz rule 

(|^) for fractional order, 

OO / 

\ ^ / —a 



where 



I''{f9) = X.[ ) f'\l'^''9). (41) 
fc=o ^ ^ 

Note that ( PD turns out to be (|36D as a tends to 1, obviously. In this respect, we 
should remark the fact that the Riemann-Liouville integrals do not define uniquely the 
Leibniz rule (0): further generalization to the rule is also possible |Q. However it is 
sufficient for the present purpose that the Riemann-Liouville integrals can fulfil (|4T|). 
With the features observed above, we accept that the Riemann-Liouville integrals yield 
the pseudo-differential operators of fractional order, 

d-" = r. (43) 

Thus we find that the pseudo-differential operator of fractional order in the Lax operator 
of the eKPi/jv hierarchy is not only an element of a non-commutative algebra but a 
concrete operator on a certain class of functions. 
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Although we do not have needed the exphcit operation of the fractional integrals 
in the derivation of the extended KP hierarchies considered in the previous section, we 
expect that the direct application of the Riemann-Liouville integral (0) is required in 
further consideration of the eKPi/7v system. For, we will need the Riemann-Liouville 
integrals when we consider the analytic solutions to the eKPi/7v hierarchies through, 
e.g., the inverse scattering method, Backlund transformation, Painleve analysis and so 
on 121 



5. Concluding remarks 

In this paper we have considered the extensions of the KP hierarchy by introducing 
fractional order integral/differential operators, the eKPi/^r- In particular, if we introduce 
the half-order integral operator, we find the resulted eKPi/2 hierarchy is analogous to 
the SKP hierarchies. Other extensions, the eKPi/jv hierarchies, are also considered and 
the fact that there exist infinite conserved currents is observed. We have also found an 
algebraic solution to the eKdVi/2 equation, the reduced eKPi/2 equation. 

To obtain the deep understanding to the eKPi/^v hierarchies, it is profitable to make 
an analysis in the Sato formulation, and also the tau-function formulation, just as in 
the case of the standard KP hierarchy and the SKP hierarchies |]T7| , |1^, |19| . Apart from 
the whole hierarchies, it will be interesting to investigate simply the integrability of 
the coupled equations like ( piop and ( pi^ ) through the Painleve analysis ||2^. Another 
approach is also attainable: the Hirota bilinear method is applicable to find special 
solutions like sohtons. As mentioned in the last section, it will be necessary to consider 
the Riemann-Liouville integrals explicitly for these analysis. 

Incidentally, apart from integrable systems there are many works on fractional order 
evolution equations for relaxation, diffusion, stochastic process and so on |^ . Although 
we have obtained the nonlinear PDE's with normal derivative through the application 
of fractional calculus, it is intriguing to formulate a systematic procedure for creating 
the PDE's with anomalous derivative, i.e., fractional order derivative. A possibility will 
be given by the application of another "Leibniz rule" in the formulation given in this 
paper, in fact, the Riemann-Liouville integrals enjoy miscellaneous types of "Leibniz 
rule" as mentioned in the last section. 



Appendix 

This appendix is devoted to present some of the important properties of the Riemann- 
Liouville integrals (p^). 

First of all, we observe the exponential law /"/'^ = 7°"'"^. For 3?ea, 3?e/3 > 0, 



1 



r(a)r(/3) J, 



dwf{w) / dz{x- zY-^{z-wf-^ 



dw{x - 
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1 

~ r(« + /?) 

= r+^/(x), (A.i) 

where we changed the integration order and used the fact that the z-integral in the 
second hne was expressed by the Beta function. This formula can be extended to the 
case when the one or both of 3fJe a and 3fte (3 is negative, by ( pSf ) and the commutativity 
of and ordinary derivative. 

Next we derive the Leibniz rule for fractional order integral/differential operator 
(PI). For SJea > 0, if we expand one of the operands, say, / in Taylor series, we find. 



1 



X 



r{f{x)g{x)) = -^ I {x-zr-^f{z)g{z)dz 
r(«) Ja 

1 r, ... i-i^^ 



Via) 



O- 1 — n 



k=0 



-1 ( 1 \ii /*x 



= E ^)f^'\^)l'^''9{x), (A.2) 

fc=0 ^ ^ 

which can be extended to 3?e a < by ordinary differentiation ( p^ ) . 

Finally we comment on the Riemann-Liouville integral of order 3ftea < 0, i.e., 
a differentiation of arbitrary order. In the definition (|38| ) the integral is well-defined, 
however if we simply put 3?ea < in the definition of /" (P7D, then the definition turns 
out to be divergent integral whatever the function / is, due to the singularity of the 
kernel [x — z)"'^^ at the upper bound. To give a well-defined meaning for the divergent 
integral, we can take the finite part of it, the Pf {partie finie) prescription. To see this 
we consider only the case — 1 < 3?ea < for simplicity, the case of 3?ea < —1 can 
be treated similarly. We now define the "regularized integral" with a positive cutoff 
parameter e as, 

I^fix) := r \x - zr-'fiz)dz 




(A.3) 



Riemann-Liouville integrals of fractional order and extended KP hierarchy 
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due to the condition /(a) = and r(a + 1) = aV{a). In the last hne of ( [A.3| ), when 
we take e — > 0, the terms of e'^^^ in the infinite sum are zero if > 1, and also the 
integral term remains finite. Thus we can perform the Pf prescription by the following 
definition, 



r(a) Ja 



- s)ds 



r{a + l) 

= r^'f'ix). (A.4) 

We therefore conclude that the definition (^) is well-defined and the Riemann-Liouville 
integral can be continued analytically to negative 3?e a. 
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